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In this paper, we investigate the properties of frustrated three-leg spin tubes under a magnetic field. 
We concentrate on two kind of geometries for these tubes, one of which is relevant for the compound 
[(CuCl2tachH)3Cl]Cl2. We combine an analytical path integral approach with a strong coupling approach, 
as well as large-scale Density Matrix Renormalization Groups (DMRG) simulations, to identify the presence of 
plateaux in the magnetization curve as a function of the value of spin S, We also investigate the issue of gapless 
non-magnetic excitations on some plateaux, dubbed chirality degrees of freedom for both tubes. 

PACS numbers: 71.10.Pm,75.60.-d 



I. INTRODUCTION 

Nowadays one-dimensional and quasi one-dimensional an- 
tiferromagnetic (AF) quantum spin systems are a very ac- 
tive theme in condensed matter physics. Thanks to efforts 
in chemical synthesis, it is now possible to obtain materials 
which can be effectively considered as one-dimensional sys- 
tems, making possible to verify the theoretical predictions. 

The natural extension of quantum spins chains are quantum 
spin ladders, which are made of two or more coupled chains. 
These ladders represent the first step between one and two- 
dimensional systems. They give rise to interesting features 1 
and have been extensively studied over the last decades. So 
far, various properties have been established both analytically 
and experimentally. 28 Going back to the problem of a sin- 
gle Heisenberg spin-S* chain, we know since the work of Hal- 
dane 9 that the chain is gapless (respectively gapped) if S is 
half-integer (respectively integer). In a similar way, there is a 
parity effect of the number of coupled half-integer spin chains 
to form the ladder. 10 A gap opens in the spectrum of antifer- 
romagnetic spin ladders with an even number of legs and their 
spin correlation functions decay exponentially. On the other 
hand, for an odd number of legs of half-integer spin, such lad- 
ders have massless excitations above their ground state and 
the decay of the spin correlations is algebraic. Experimen- 
tal investigations have confirmed these predictions, for exam- 
ple the observation of a gap in the spin-1/2 two-leg ladders 
Si'Cu20 3 n andCu 2 (C 5 H 12 N 2 ) 2 Cl4. 12 

New properties arise when the role of transverse boudary 
counditions is taken into account. The results quoted above 
are valid only for ladders, which correspond to open bound- 
ary condition (OBC). Applying periodic boundary conditions 
(PBC) in the rung direction to form a spin tube seems to cause 
the opening of a gap in both even and odd cases. For an even 
number of legs, this is explained in terms of the formation of 
spin spinglets in the transverse direction. The reason is dif- 
ferent in the odd case, where the PBC induce geometrically 
frustrated interactions. The direct consequence of this frustra- 
tion is a twofold degenerate dimerized ground state with an 
excitation gap above it. 13 

Among the questions emerging from the study of these 
quasi one-dimensional systems, an important one concerns 



their magnetization process when an external magnetic field 
is turned on. Classicaly, the magnetization curve of such sys- 
tems is expected to be a straight line until the saturation field. 
But at low enough temperatures, quantum effects begin to play 
a role and magnetization plateaux can appear. This has been 
observed in various chains and ladders spin systems. 14 ~ 17 

A condition, neither sufficient nor necessary, for the exis- 
tence of magnetization plateaux in a quantum spin-S* chain 
has been found by Oshikawa, Yamanaka and Affleck. 18 This 
condition, which was later extended to ladder systems, 6 re- 
stricts the possible values of the magnetization for a plateau. 
It reads 



N(S-m)e 



(1) 



where m is the on-site magnetization and N the number of 
spins per unit cell. This result has been obtained through 
a generalization of the Lieb-Shultz-Mattis (LSM) theorem. 19 
Although the low-energy state constructed in this approach 
has the same total magnetization as the ground state, they 
present bosonization arguments indicating that, in general, 
low-energy states appear also in different magnetization sec- 
tor. In 2009, Tanaka, Totsuka and Hu (TTH) have used a spin 
coherent-states path integral approach to recover this condi- 
tion. 20 The main advantage of their method is that it can be 
applied for any value S. 

In this paper we investigate successively the effect of a 
magnetic field on two different types of three-leg spin tube 
of spin-S*, namely the simple spin tube and the twisted spin 
tube, which will be describe below. For each one, we com- 
bine analytical and numerical methods and proceed as fol- 
lows. First we apply to the spin tube the TTH approach, which 
take into account the effects of the Berry phase appearing in 
the partition function. This leads to an effective field theory 
for the spin tube, from which a condition on the magnetiza- 
tion plateaux is infered and we give an estimate the region of 
existence of the plateaux. We also discuss the question of hav- 
ing gapless non-magnetic excitations but also several different 
gapped phases for these excitations. Then, in the case of half- 
integer spins we study the limit of strongly coupled chains. A 
new non-magnetic degree of freedom appears in this regime 
for the magnetizations of the lowest and highest magnetiza- 
tion plateaux, namely a right or left chirality. It comes from 



the twofold degeneracy of the ground state. In this limit we 
show the possibility of the existence on those plateaux of a 
quantum phase transition. As the coupling along the chains 
is increased, the chirality degree of freedom may go from a 
critical to a gapped regime. This behaviour has recently been 
observed in the simple spin tube of spin- 1/2 by Okunishi et al. 
using density-matrix renormalization-group (DMRG) calcu- 
lations. 21 We expect this transition to be well accounted by a 
pertubated XXZ effective Hamiltonian and support this state- 
ment by a very simple qualitative numerical result. Then we 
perform DMRG calulations to study the magnetization pro- 
cess. Analyzing the entanglement entropy and the local mag- 
netizations, we finally examine the different phases occuring 
for the chirality on the plateaux. 

The rest of the paper is organized as follows. In Sec. II, 
we consider the case of the simple three-leg spin tube by us- 
ing the path-integral approach to understand the appearance 
of some magnetization plateaux, then we investigate the role 
of the chirality degrees of freedom, and finally we compare 
our predictions to DMRG simulations in the 5 = 3/2 case. 
In Sec. Ill, the case of the twisted spin tube is addressed fol- 
lowing the same strategy. Finally, we draw some conclusions 
and discuss possible perspectives in Sec. IV. Some technical 
details about the duality transformation are given in appendix. 



II. SIMPLE SPIN TUBE 

A. The model 

The Hamiltonian of the simple spin tube, which is the first 
geometry that we consider, reads 



H = H x + #|| + H h 

j a=l,2,3 

-^ii = ^ii Z^ Z^ s a j.s a j + i 

j <*=1,2,3 



(2) 



j a=l,2,3 

where S a ,j is the spin-5 operator, Jn > is the intrachain 
AF coupling, J± > the AF rung coupling and h the mag- 
netic field along the z axis (Fig. 1). The subscript i (a) repre- 
sents the site number in the chain (rung) direction. The tube 
structure (PBC in the rung direction) induces frustration in 
this simple nearest-neighbour model. Applied to this model, 
the OYA condition ( 1 ) takes the form 



3(5 -m) eZ, 



(3) 



as there are three spins per unit cell. 

So far, the Hamiltonian (2) has already been investigated 
in previous works. By using bosonization techniques for the 
spin 5 = 1/2 and in zero magnetic field, Schulz has suggested 
that the tube has a spin gap induced by the geometric frustra- 
tion. 22 DMRG calculations have confirmed this statement by 




FIG. 1 . Lattice structure of the simple spin tube. 



establishing the existence of a spin gap, 13 which is explained 
in terms of a two-fold degenerate dimerized ground-state with 
broken translation symmetry. New questions arise when one 
of the couplings between the rungs is changed, making the 
tube asymmetric. In that case, it seems that the spin gap van- 
ishes for a small but finite asymmetry. 23 The effect of a mag- 
netic field has also be addressed. 24_28 For higher half-integer 
spins, only a few results are available. In the 5 = 3/2 case, 
the existence of a spin gap, as for 5 = 1/2, has been reported 
recently. 29 Moving to integer spins 5, the simple spin tube 
displays interesting properties. For weakly coupled chains, 
another parity effect has been established for the low-energy 
properties using the Non-Linear Sigma Model (NLSM). For 
an odd number of legs, the lowest magnon band of the model 
(2) is six-fold degenerate, compared to the three-fold degen- 
eracy of non-frustrated systems, namely tubes with an even 
number of legs or ladders. 3032 Turning on a uniform or non- 
uniform magnetic field, Sato predicted a succession of quan- 
tum phase transitions between critical phases as the field is 
increased, along with again an even-odd effect. 30 Introduc- 
ing an asymmetry in the rung couplings, a NLSM analysis 
and DMRG results have shown the existence of 25 quantum 
phase transitions between gapped phases when varying the 
anisotropy parameter. 33 

Finally, we discuss the classical configurations of the spin 
tube (2). For decoupled triangles with no external field, the 
three spins are simply lying in the same plane with angles of 
27r/3 between them. If the triangles are now coupled by the 
J|| term, the situation does not change as the longitudinal cou- 
pling is not frustrating and satisfied with a fc|| = n order in 
this direction. A magnetic field arrange the spins in an "um- 
brella" configuration, where the three spins on the triangle are 
equally polarized and have angles of 27r/3 between their pro- 
jections in the plane. Thus, the classical ground state of the 
simple spin tube is simply an umbrella configuration on each 
triangle with a canted order along the tube. We parametrize it 
as 



£>a, 



'(-l)->sm(# )cos(<^ 
S | (-l)i s in(0 o )sin(^ 



(4) 



cos(6>o) 



where cos( 



= sojx-hjii) and & = (« " 1)2tt/3. 
This state breaks the U(l) symmetry around the z axis, one 
site translations and parity transformations. More precisely, 
the latter one is a symmetry which is going to play a very 
important role in this system. It is related to what we dub 
the chirality degree of freedom. Consider the chirality vector 



order parameter: 



Xj ~ q /.("a,j x ba+l,j) , 
6 a =l 



(5) 



which is invariant under cyclic permutation of the three chain 
indices (i.e. translations in the transverse directions) but 
changes sign under the permutation of two chains. For the 
classical configuration (4), it reads \j &■ Ea s ' n ((rf — 
Pa+i)/ty- Thus, the choice of a given classical configuration, 
namely the choice of (p^ — (p^+i — ±2-7r/3, fixes the sign of 
this order parameter and breaks the Z2 symmetry explicitly. 
This in turn will have important consequences in the analysis 
of the chirality behaviour within the path integral approach. 



B. Path integral approach 

1. Derivation of a low-energy action 



is the coherent-state expectation value of the Hamiltonian, or 
simply the Hamiltonian for the classical configuration. The 
other term is the Berry phase one and corresponds to the sur- 
face area (or the solid angle), /dr[l — cos(#(t))]9 t <p(t) 
in spherical coordinates, enclosed by each spin during its 
imaginary-time r evolution. We want to build a low-energy 
effective theory from this starting point. The method consists 
in finding the classical ground state of the system and then 
adding the quantum fluctuations to derive an effective action. 
We start from the classical ground state discussed in 
Sec. II A and now we add the fluctuations around this state, 
writing 



Va -> V 



<x,3 




0n + SO. 



fa,. 



01,0 
(a- 



2tt 



(6) 



<Po 



and expand the spin components up to second order in 58. The 
calculation of the SU(2) commutation relations [S^^Sq] 
leads to introducing a new set of variables Tl a j, defined by 



We begin the study of the simple spin tube (2) following 
the method recently developped by Tanaka, Totsuka and Hu. 20 
They used a Haldane's spin coherent-state 34 path integral ap- 
proach to re-derive the OYA condition (1) for the Heisenberg 
chain with easy-plane single-ion anisotropy. The interest of 
the method is its validity for any value of the spin S. Hal- 
dane's analysis leads to an action comprising two terms. One 



n„, 



-s 



sin(e )69 aJ + -cos(9 )Se 2 a ^ 



(7) 



which are the conjugates of the (p a j's. It ensures to have the 
correct commutators for the spin operators. Then we rewrite 
these operators as functions of the conjugate fluctuations vari- 
ables ip a j and n aj as 



S Zj ~ (-I) i e ±i [ (a ~ 1) ^ +V '-']S 



s z a , « scos(0 o ) + n c 



' . , ft s rn 1 S 2 _l 2 

Smm S 2 sm(6 ) a ' j 2 S 2 - m 2 S^m(9 ) a ' j 



(8) 



where m — Scos(6o) is the classical magnetization per site. 
Casting these expressions into the action, taking the contin- 



uum limit and keeping terms up to second order in the fields, 
we obtain the low-energy effective action 



S[{U a },{<p a }] = j ' drdxl Yl 



= 1,2,3 

1 m 2 



\a,]\\{S 2 - m 2 )(d xVa f + a Uj^ + lj x -^— - 



n 2 



J± S 2 - m 2 



^-M 1 ' W 2 (nina + nana + nanx) 

\ 2 b z — m z J 4 a 

/q 

—mj ± [u.i(ip 3 - tp 2 ) + n 2 (yji - (£3) + n 3 (^ 2 - <Pi)} + i ^ 



[(<fl ~ V2) 2 + O2 - <yS3) 2 + (<f3 ~ Pi) 2 
S — m s 



a=l,2,3 



d T (p a - U a d T ip 



(9) 



with a the lattice constant. We see that all the fluctuations, 
transerve or longitudinal, are coupled. The last two imagi- 
nary terms come from the Berry phase part of the action. It 
is important to stress that the d T ip a terms, although being to- 
tal derivatives, can not be dropped. Indeed, the fields tp a are 



angular variables defined on a circle and thus this term counts 
the winding number of each field. 

As at this order the action is gaussian in the fields H a , we 
can integrate them out and the action becomes 



£[{&*}] = Sch[4>i,4>2\ + S a [(j> a ] 
Sch[<h,<h] = J dTdxl^ 2) [(drfa) 2 + (drfo) 2 ] + \\^ 2) [(d x <f>i) 2 + (d x <f> 2 ) 2 ] +m 2 (^ + , 

S 8 l<t> s ] = I drdx^X^idr^) 2 + ^{d^sf+iZ^-^d^ 



V 



>a) - in((f>id T <f>2 



</>2<9r0l) 



(10) 



where S c h denotes, for reasons which will become clear later, 
the chirality part of the action and S s the symmetric one. We 
have made the change of variable <p = U0, where 



U = - 



and have rescaled the symmetric field as 
coefficients of the action (10) read 





A^ 1 . 2 ) = 

a UJ|| - 

\£>V = aJ l{ (S 2 

M 2 =3JiiJ_L- 



1 



3/ ^ 
2 JJ - S 2 



, 4« 



s) 



a(4Ji| +3J ± )' 



'■), A^=3aJ||(5 2 -m 2 ), 



a 4 J,, 



3 t m 2 

2 J - L S 2 -m' 



W 



a 4 Jii 



fJ± 



S 2 -m s 



(12) 



The symmetric field is now decoupled from <f>\ and <p2, and 
we will study them separately. 



2. Symmetric action and magnetization plateaux 

We notice that the action S s for the (f> s field obtained in (10) 
has the same form than the action of the Heisenberg chain 
in a magnetic field. 20 The term id T (f> s comes directly from 
the Berry phase part of the action discussed above. In order 
to understand its role on the low-energy physics, we apply 
a duality transformation 35 on this action (details are given in 
Appendix A). The dual action finally reads 

1 



5[*J= / drdx{ ^ S (V<I S ) 2 



gicos I 2-7T 



- „ 5 — m 

$ s + 3 x 



(13) 



where $ s is the dual field, K s = l/\J \ { T s) \ x s \ V = (d T , d x ) 
and g\ is a constant we have not computed. 

We find a gaussian action perturbed by a cosine term. It 
is the sine-Gordon action in (1+1) dimension plus a spatial 
modulation 27r3(S'— m)x/a of the cosine. This modulation is 



a direct consequence of the topological Berry phase. We can 
therefore separate two cases. If 3(5 — m) ^ Z, the cosine is 
incommensurate and will average to zero by integrating over 
space. It remains a simple gapless gaussian model and there 
will not be any plateau in the magnetization curve for general 
values of to. We can understand the effect of the Berry phase 
in terms of protecting the system from the vortices. On the 
other hand, if 



3(5 -m) € 



(14) 



the cosine is commensurate and we recover the sine-Gordon 
model. A gap can open in the spectra, causing the emer- 
gence of plateaux in the magnetization curve for these par- 
ticular magnetizations. For example, we expect to observe a 
plateau for average magnetization per site 1/6, 1/2, 5/6 and 
7/6 in the S = 3/2 case. 

More precisely, the presence of the cosine term is not 
sufficient to open a gap. It has to be relevant in the 
renormalization-group (RG) sense and it will depend on the 
microscopic parameters. From the well-known action (13) 
(without modulation) we compute the scaling dimension A 



A = 7T 



\ 



3Jji(5 2 -to 2 ) 



J± 1 



Mi 

3J, 



(15) 



For the cosine to be relevant A has to be smaller than 2 
and thus we expect plateaux in the weakly coupled triangles 
regime. 

We end this discussion with some comments. First, it 
is worth noting that the condition (14) does not predict the 
to = plateau found numerically by DMRG. 13 But in the 
action ( 1 3), higher harmonics of the cosine term have been 
dropped, especially the first one which actually predicts this 
plateau. However, as the harmonics terms would be less rel- 
evant than the fundamental, we can not conclude about the 
to = plateau with this analysis. In fact, the zero-field 
problem needs a different approach taking into account the 
0(3) symmetry of the model. 33 Second, the action (13) is the 
same as the one obtained for spin-1/2 using bosonization 18 , 
but coming here from a large S method. So, with the condi- 
tion (14) we have simply recovered the OYA result (3) on the 
possible values for magnetization plateaux to occur. 

Finally, it is important to note that the factor of 3 in front 
of the Berry term s ~ m d T (f) s is not an artefact of the transfor- 
mation (11) and of the rescaling of the field. The rescaling is 
necessary for the following reason. As the (^fields are angular 
variable, they satisfy 



(fi(T = p) = <pi(r = 0) + 2nn 2 



(16) 
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FIG. 2. (Color online) Scaling dimension of the cosine operator in 
the action (13) as a function of J\\/J± for S = 3/2 at the values of m 
fulfilling the condition (14). The perturbation is relevant and opens a 
gap in the spectrum if A < 2, thus the plateaux should disappear for 
J\\/J±> 0.1. 



and the same condition for spatial periodic boundary condi- 
tions in x = and x = L. That gives similar conditions 
for the new fields <fi, for instance the antisymmetric combina- 
tion </>i satisfies 4>i(f5) = </>i(0) + 27r(rt2 — ni)/\2. But, 
given the form of the action S c h for <j>\ and </> 2 (see below), 
we expect those fields to be small and therefore to have no 
winding, i.e n\ = n 2 = 0. It remains for <f> s the periodicity 
4> S {P) — </> s (0) + 27rri3/v3- Thus the rescaling <p s — >• <p s /\/3 
is required to ensure a correct 27r-periodicity. 



3. Chirality degree of freedom 

We now focus on the action S c h for the fields <\>\ and 02, 
which stands for the chirality (non-magnetic) degree of free- 
dom we will introduce in Sec II C 1. Before going into the 
technical details, a first comment is in order here. As we 
mentioned before, the chosen classical configuration on top 
of which the path integral approach is constructed explicitly 
breaks the Z 2 symmetry related to chirality. However, from 
the arguments shown below, we expect it to be able to repro- 
duce almost all phases potentially observable for the spin tube. 

It is convenient to introduce the two complex conjugate 
fields f = ^i+ i<t>2 and W* — <j>\ — icj>2 and, after rescal- 
ing the time, we get the action 



5[*,** 



drdx 



xiwP + m 2 !*! 2 



2/z(**d T *-*d r **) 



where K = \ X ( T ' ' \ { x ' ' and M 2 = M 2 



(17) 



We write 



expansion and higher order terms are expected in the general 
action. Although there is a mass term, we believe that the 
Berry term (j,(if?* d^ — *f?d T if?*) may induce strong effects on 
the behaviour of ty, namely the possibility to have a gapless 
phase. We propose to treat qualitatively this question by using 
the symmetries to write a general action including important 
higher order terms. 

Going back to the initial fluctuation variables defined in (6), 



it is instructive to rewrite the field if? as <£ 



r(ip 3 



uj 2 tf2), where u> 



So under a circular permutation 



of the sites on a triangle, which let the system invariant, the 
field ^ grabs a phase factor if> -> cjif?. Then the most general 
action invariant under such transformations has the form 

S[f, **] = / drdxl -iv~|W| 2 - /i(**a r * - *<9 r **) 



+ /3* 3 + /r** 3 + /(|*u s ) + 



(18) 



Writing ^> = pe 10 , the only potentially gapless degree of free- 
dom is the phase field 9, and the most general action reads 



S[6] 



f drdxl ^K\V9\ 2 



A 3 cos(36>) + A 6 cos(66>) 



+ i^d T 9 



(19) 



(...) to remind that we are currently working in a second order 



where K, A3, \§ and fi are phenomenological parameters. 
The important point to note is that the original Berry phase 
causes the i/j,d T 9 term that forbid the vorticity (see the discus- 
sion for the symmetric action). 

The effective action (19) tells us that there are four possible 
phases : 

i) When the stiffness K of the field 9 is large enough, the 
scaling dimensions A3. 6 ~ 1/ K of the cosines are small and 
they are relevant so that (if?) ^ 0. We first assume that A6 has 
the same sign as A3. In this case we have only three equiva- 
lent solutions for ("f?) in which two of the three fields ipi have 
the same value. To understand the consequences of this, let 
us go back to eq. (9) and rewrite the before the last term as 

ni(^3-^2)+n 2 (v3i-^3)+n 3 (</?2-^i) oci(n**-n**), 

where the complex conjugate fields II and II* have the same 
definition than the fields ^, ^?* but with respect to the fluctua- 
tions II Q . Having (W) 7^ implies (II) 7^ 0, which translates 
into a homogeneous renormalization of the classical value of 
the magnetization. This correction is equal for two spins but 
different for the third one (note that the total magnetization 
is kept unchanged). We call this phase the symmetric spin 
imbalance. 

ii) The cosine operators are again relevant so that (\P) 7^ 
but now A6 has the opposite sign than A3. As for the phase 
(/) there is an homogeneous spin imbalance but with three dif- 
ferent values for the magnetization, and we dub this phase 
an asymmetric spin imbalance phase. The transition from the 
threefold degeneracy of the phase (/) to the sixfold degeneracy 
of this phase corresponds to the double sine-Gordon model, so 
it belongs to the universality class of the Ising transition. 36 



Hi) For a sufficiently small stiffness, the cosine opera- 
tors become irrelevant. In that case we perform a duality 
transformation as previously to take into account the role of 
the \d T 9 term. We end with two terms (1/A")(V9) 2 and 
A2COs(27r(8 + rjx)), where 9 is the dual field of 9 and 77 a 
phase modulation, a priori function of the microscopic param- 
eters. For general 77 phases, the cosine is not commensurate 
and that eventually leads to a gaussian model. This phase is 
characterized by ( , 5) = and algebraic correlation functions 
for the 9 field. We have a conformal field theory with central 
charge (see Sec. II D 2) c = 1. 

iv) (\fr) = and the correlation functions for the 9 field 
are short -ranged. Here some comments are in order concern- 
ing the action (19). The last term, which originates from the 
Berry phase, has the effect to suppress vortex configurations 
for the field 9 and is responsible for the gapless phase (Hi). 
It is also the same scenario found to explain that the field <p s 
is gapless in general except for particular values of the mag- 
netization. In the present case, as long as the parameter ji 
has a generic value, we said that the cosine operator of the 
dual field is forbidden. 37 Thus the transition from the phase 
(Hi) to this short-ranged phase should not be in principle via a 
Berezinsky-Kosterlitz-Thouless (BKT) transition. It would be 
rather because K -> 0, much in the same manner than a XXZ 
chain enters into the ferromagnetic phase when the exchange 
anisotropy parameter becomes sufficiently negative. Such a 
ferrochiral phase has been found for S = 1/2 and weakly 
coupled chains in a wide range of the magnetic field by Sato 
in Ref. 27. He argued that this order should also survive when 
entering the plateau state for a moderately larger rung cou- 
pling. 



doublets. These states are 
1 1 L) = -L (| m) 

\i L) = -^ (| Ut) 



Ut) 



in)) 
m» 



(20) 



1 1 R) = -7= (I m) + w _1 | ut) + "I in)) 

UB> = i(l ut) + ur 1 ! m) + u\ m» , 

where u — e 1 ^ . The indices L and R represent the chirality 
and f, 4- the z-axis projection of the total spin of the triangle. 
It is important here to make the link with the field W oc (^3 + 
uitpi + uj 2 (p2 defined in the path integral approach. Given the 
similar form of the states (20), it indicates clearly that the field 
^f (or equivalently the fields <\>\ and (j>2 ) describes this chirality 
degree of freedom. This is also consistent with the fact that 
the Berry phase of S c h m (10) disappears in the opposite limit 
Jl/J|| -»• 0. 

In this strong coupling limit, we keep only the four states 
(20) to describe the low-energy physics around the zero mag- 
netic field level crossing. To first-order in J11 / J± the effective 
Hamiltonian reads 

3 3 

(21) 
where Sj is the triangle total spin-1/2 operator. We define 
the pseudo-spin- 1/2 chirality operators r . They exchange 
chiralities L and R such as 



■L)=0, 



t-\-L) = \-R), 



t+\-R) = \-L), t-\-R)=0. 



(22) 



C. Strongly coupled chains : effective models and chirality 



1. First-order perturbation Hamiltonians 



In Sec. II B, we have discussed the model (2) regardless 
of the value of the spin S or the strength of the coupling pa- 
rameters J11 and Jj_. We eventually found a condition on the 
magnetization plateau values. In this section, we focus on 
the half-integer spin case and on the strong coupling limit be- 
tween the chains, or weakly coupled triangles, J±/J\\ —t 00. 
For a given half-integer S, we study the lowest and highest 
magnetization plateaux, namely m = 1/6 and m = S — 1/3, 
as we will show they both can be described in terms of an 
additional chirality degree of freedom. 

We start with the S = 1/2 case in the strong coupling limit, 
and consider first the extreme case Jji = where the system 
is made of decoupled triangles. The ground-state of a trian- 
gle is fourfold degenerate at h = 0, with two chiral spin-1/2 



By construction, this effective Hamiltonian describes the sys- 
tem from zero magnetization up to the first plateau m = 1/6 
with (Sf) = +1/2, where only the two polarized states re- 
main. This model has been studied both analytically and nu- 
merically. 13 ' 22,25 Its spectrum displays a small plateau at mag- 
netization m = 0, the spin gap arising from the dimerization 
of the ground state as explained in Sec. I, A strong enough 
magnetic field closes the gap and the system is then described 
by a two-component Luttinger liquid, with both spin and chi- 
rality modes being gapless. 25 Increasing again the magnetic 
field drives the system to the magnetization plateaux where 
only the two S z = +1/2 states are present. The chirality is 
described by the XY Hamiltonian 



tfe//=^£[l + 4(r; 



'j+i 



T 3 T t+ 



i)]- (23) 



and is then critical. 

This description remains valid for higher half-integer spins 
S. The low-energy space of one triangle at zero magnetic field 
is always spanned by two degenerate chiral doublets whose 
spin projections are S z = ±1/2, so the above description 
used to derive the effective Hamiltonian can be repeated. The 



region from zero field up to the first plateau m = 1/6 is 
described by an Hamiltonian of the form (23), with only a 
change in the numerical constant of the chirality operators. 
On the plateau, the physical spin is frozen to +1/2 and the 
chirality is governed by an XY model which reads 



H, 



' if ~ 12 5J 

3 



a{ T] Tj+1 



'j ']+ 



i)]> (24) 



where the single parameter is a = (25 + l) 2 . 

Starting again from the decoupled case Jn = 0, we ob- 
serve that the above chirality description can also be used for 
the highest magnetization plateau m = S — 1/3, where the 
isolated triangle ground state is also twofold degenerate. Us- 
ing again first-order perturbation theory, we find the chirality 
states on the plateau are also given in terms of a XY model of 
the same form as (24). 



2. Range-2 CORE Hamiltonians 

However, it turns out that these first-order effective Hamil- 
tonians do not capture entirely the behaviour of the chirality 
on these two extreme plateaux, as we observed numerically by 
measuring the central charges the existence of gapped phases 
for some range of the coupling (see Sec. II D 2). 

A way to go beyond the simple first-order perturbation 
theory is to use a Contractor Renormalization 38 (CORE) ap- 
proach to compute numerically an effective Hamiltonian. The 
CORE technique is a non-perturbative method of renormal- 
ization in real space for lattice systems, used to build effective 
Hamiltonians reproducing the low-energy physics. It has been 
shown to give quantitative results for instance for various an- 
tiferromagnetic models 39 , including the presence of magnetic 
field. 40 Here, we truncate the calculation of the effective in- 
teractions to range-2, i.e. we consider only two coupled trian- 
gles. This is a quite simple computation, but it already gives a 
qualitative improvement over lowest-order perturbation, since 
it produces an effective Hamiltonian of the XXZ type : 
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eff 
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(25) 



We show in Fig. 3 the computed values of the parame- 
ters J xy , J z and their ratio A = J z /J xy governing the be- 
haviour of the model (25), for the plateau rn = 7/6 in the 
case S = 3/2. For coupling values J\\/J± < 0.042, the 
system is in the regime | A| < 1, the gapless XY phase with a 
central charge c = 1. As the coupling is increased, we see that 
the negative J z component decreases and at the critical value 
J\\ C /J± = 0.042, the system enters in the regime A < — 1, 
corresponding to the ferromagnetic ("ferrochiral" here) phase 
c = where all the triangles have the same chirality L or R. 
Note that the change of sign of J xy does not change the nature 
of the phase as only the sign of J z , which remains negative, 
is important. Also, we have checked by exact diagonaliza- 
tion (ED) for larger system lengths (up to L = 10) that there 



is a level crossing close to this critical coupling. The quan- 
tum numbers of the ground state are compatible with the XY 
to ferrochiral scenario. Finally, the same scenario for A oc- 
curs for the lowest plateau m = 1/6 at a critical coupling 
J\\ C /J± = 0.256 (however we will see with the DMRG re- 
sults that for this value we are no longer in the plateau phase). 
Even if this approach is straightforward, it allows to explain 
qualitatively the possibility of a phase transition from a critical 
to a gapped phase. On the quantitative side, although longer 
range effective interactions are expected to play a role (see 
below), the critical value J\\/J± < 0.042 found here is very 
close to the transition value observed with the DRMG. 




FIG. 3. (Color online) Values of the effective XXZ Hamiltonian pa- 
rameters J xy , J z and A = J z /J xy describing the chirality on the 
m = 7/6 plateau for a spin tube S = 3/2. The inset shows that at 
a critical coupling J\\/J± = 0.042 the chirality excitations become 
gapped and the model enters the ferrochiral phase. 

As this non-perturbative CORE computation only involves 
solving two triangles, we can also treat higher spins. The Ta- 
ble I shows, for different values of the half-integer spin S, the 
critical values jf, c and J,j l c at which the chirality on the low- 
est and the highest plateaux is expected to undergo a quantum 
phase transition from the XY phase to the Ising phase. For 
both plateaux, we observe the gapless phase shrinks in the 
large S limit. Again, the results given by this simple method 
are encouraging and we believe them to be actually quite accu- 
rate. Indeed, in the spin-1/2 case we find the value 0.500, very 
close to the value 0.496 computed with DMRG where the chi- 
rality enters the ferrochiral phase 21 . In the following we will 
present DMRG data supporting these results for S = 3/2. 





5 = 1/2 


5 = 3/2 


5 = 5/2 


5 = 7/2 


5 = 9/2 


r l 


0.500 


0.256 


0.157 


0.108 


0.079 


jh 

J \\,c 


0.500 


0.0420 


0.0140 


0.0071 


0.0042 



TABLE I. Critical values J« c and Jii l c of the transition from an XY 
effective model to a ferromagnetic Ising one to describe the chirality 
behaivour on the lowest and the highest magnetization plateaux. 



3. General effective Hamiltonians 

Recently, Okunishi et al. have derived the second-order ef- 
fective Hamiltonian in the spin-1/2 case 21 , which is a special 
case since it has only one plateau. New terms appear, such as a 
negative t z t z one which can drive the system into an ordered 
phase. This is in agreement with the CORE calculation. Our 
goal is to propose an effective Hamiltonian capturing these 
phases and the transition. As the argument of the phase factor 
u> in Eq. (20) is nothing else than the transerve momentum, the 
effect of the operators r is simply to shift the triangle mo- 
mentum of ±2-7r/3. By transverse momentum conservation, 
the most general effective Hamiltonian, on both the first and 
last plateau, can be written as 



H eff = 
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2 V 
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J Z Z 
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a transition to the ferrochiral phase when the magnitude of 
the negative A becomes sufficiently large. We associate this 
second phase to the case (iv) predicted by the path integral ap- 
proach. The first case eventually leads to the appearance of 
a gap in the chirality degrees of freedom in favor of a spin 
imbalance phase similar to the one found in Ref. 21 but with 
one-step breaking of the translation symmetry. 

Indeed, having a non-zero expectation value for the field x) 
make that the operator t x ~ (— l) x cos(x) is also non-zero. 
This correponds to a symmetry breaking, as this operator is 
directly related to an imbalance in terms of the original spin 
operators. 42 As this phase is incompatible with a non-zero av- 
erage magnetization in the z direction, it is inaccessible to the 
path integral approach we have presented before. Also, the 
critical point at which the spin imbalance phase would oc- 
cur is expected to be very close to the regime in which the 
ferrochiral phase appears (which, in the absence of the three- 
body t + t + t + + h.c. term, occurs for A = —1). We the 
believe this phase to be only present in a very narrow range 
between the XY and the ferrochiral phases. 



(26) 



where we have dropped, for example, second-neighbour ex- 
change terms. The values of the parameters of this model have 
been calculated up to second-order in J\\/J± (see 21). The 
T + T + T + term was ]3 V i 0US iy absent in our range-2 CORE cal- 
culation but we would except it to appear for a higher range 
one. Notice that in this language the Z2 symmetry associated 
with chirality is just r z — > — t z which can be obtained for ex- 
ample with a rotation of it around the x axis. At this point one 
can try to make connection with the results of the path integral. 
We have to keep in mind that, by construction, the Z2 chiral- 
ity symmetry is broken within the path integral approach. This 
corresponds to placing the effective spin chain above at a non- 
zero average homogeneous magnetization (t z ) ^ 0. 

The bosonized form of this effective Hamiltonian was writ- 
ten recently in Ref. 2 1 as 



dx 



1 

2k 



(d xX ) 2 + 2n{d x xf 



H = 



+Ai cos(2v27rx) + A2 cos(6v27rx) 



(27) 



where v is a Fermi velocity, n is the Luttinger parameter and \ 
is dual to x- The second cosine operator is radiatively gener- 
ated by the presence of the t + t + t + + h.c. term. The first co- 
sine operator is relevant for k < 1 while the second becomes 
relevant for k > 9. The gapless phase obtained between these 
two critical points is associated with the case (Hi) predicted 
by the path integral in Sec. II B 3. For J\\/J± —> 0, the effec- 
tive Hamiltonian (26) reduces to an XY model corresponding 
to K = 2, thus both cosine terms are irrelevant and the chi- 
rality is in a critical phase. When the coupling ratio J\\/J± 
is increased, it turns out that k increases too as we are in the 
ferromagnetic regime, so the first cosine cos(2v27rx) will al- 
ways be irrelevant. 41 We predict the gapless phase to disap- 
pear and two gapped phases should appear successively. The 
first one is caused by the cos(6-\/27rx) term becoming relevant 
for some critical negative A c , and the second corresponds to 



D. DMRG results for S = 3/2 

1. Magnetization plateaux 

In order to verify the previous predictions about the mag- 
netization plateaux, which were established thanks to large-S" 
techniques, we have performed numerical simulations using 
the DMRG algorithm 43 for S = 3/2 spin tubes with open 
boundary conditions (OBC) along the legs. We consider sys- 
tem lengths up to L = 64. Typically, we kept up to 2000 states 
and perform 20 sweeps, which is sufficient to have a discarded 
weight smaller than 10~ 8 or less. 

In Fig. 4, we plot a typical magnetization curve obtained in 
the strongly-coupled chains regime for J\\/J± = 0.1. Large 
plateaux are observed below saturation for magnetization per 
site to = 1/6, 3/6, 5/6 and 7/6, which correspond to the 
condition (14) that we have found with the field theory. Re- 
member also that an to = plateau was predicted, but it has a 
different nature (dimerization of the ground-state), and on the 
scale of the figure, it is hardly visible. 

In order to map out the phase diagram, we perform a finite- 
size analysis of the widths of each plateau for several cou- 
plings. Resulting data are shown in Fig. 5. While it confirms 
that all the plateaux found for J\\/J± = 0.1 are present in 
the thermodynamic limit, we do observe that each of them 
disappears for some different critical ratios of the coupling 
constants. These critical values of J\\/J± are summarized in 
the Table II along with the predicted values coming from the 
formula (15). It is important to mention that because the tran- 
sitions between the plateau phases and the gapless phases are 
expected to be of the BKT type, it is difficult to locate them 
accurately. We also note that the predicted values are in a 
roughly good agreement with those coming from the DMRG. 
The main qualitative difference is that the path integral ap- 
proach predicts that the highest plateau, m = 7/6 here, should 
be the most robust. This prediction of the plateaux disap- 




FIG. 4. Magnetization per spin curve (ra) vs magnetic field h for 
the simple spin tube in the case S = 3/2. DMRG simulations were 
performed with L — 64 and Jii/Jj. = 0.1. Finite-size steps are 
almost not visible on this scale. 
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FIG. 5. (Color online) Finite-size scaling of the plateaux widths for 
different couplings Jv, /J±. 



pearing gradually with to seems to be a general feature of 
the path integral approach for chains or coupled chains. By 
performing extensive simulations, we arrive at the phase di- 
agram shown in Fig. 6, where we indicate both the plateaux 
phases and the chirality phases discussed in the next subsec- 
tion. It is clear that the predicted disappearance pattern for the 
plateaux is not recovered exactly, but we have to remember 
that the result (15) comes from a large-5 approach. We plot 
the saturation field h sat , which can be found analytically to 
be hsat = (3Jj_ + 4J||)S r . We also indicate the existence of 
the to = plateau, which we observe on the whole range of 
J\\ J J± but is not visible on this scale. For instance, we find 
the spin gap to be of order 2.10~ 2 J± for J\\/J± = 0.1, a 
value in good agreement with other DMRG calculations. 29 





m= 1/6 


m = 3/6 


m = 5/6 


m = 7/6 


rDMRG 
Hi c 


0.20-0.22 


0.14-0.16 


0.16-0.18 


0.22-0.26 


1 PI 


0,066 


0.074 


0.098 


0.19 



TABLE II. Critical values J« c for the plateaux, in unit of J±, for 
the spin tube with S — 3/2. We indicate both the values found by 
DMRG and with the path integral. 




FIG. 6. (Color online) Phase diagram of the three-leg spin tube with 
S = 3/2 as a function of the coupling Jm and magnetic field h 
(Jx = 1 is the unit of energy). Several magnetization plateaux 
can be observed (filled areas) and an additional m = plateau 
is found (bold line). Data correspond to numerical simulations on 
3 x 32 lattice with DMRG. Inside the extreme plateaux m = 1/6 
and m = 7/6, hashed red regions correspond to critical chirality 
phases (see Sec. II D 2). 



2. Entanglement entropies and conformal charges 

Since the chirality degree of freedom has been shown to 
emerge for the extreme plateaux for half-integer spin, we 
now use large-scale DMRG simulations to investigate it. For 
S = 1/2, chirality is only expected on m = 1/3 plateau, and 
has already been confirmed numerically. 21 In our paper, we 
consider the next case, i.e. S = 3/2. 

In order to check the existence of a chirality phase transi- 
tion on the extreme plateaux, we simply compute the block 
von Neumann entropy S v n (£) which exhibit two different be- 
haviours for large blocks £ and OBC: S v n(£) saturates to a 
constant when the system is fully gapped, whereas S v n (£) — 
(c/6)log£ where c is the central charge of the underlying 
conformal field theory. 44 In order to minimize finite-size ef- 
fects, we will consider the conformal block length d(£\L) = 
(L/tt) sin(£n/L). 

Guided by our CORE analysis, we first choose a small cou- 
pling J|| / Jj_ = 0.02 where chirality is expected to be gapless 
on both m = 1/6 and to = 7/6 plateaux. As it is shown 
in Fig. 7(a), numerical data are compatible with a gapless be- 
haviour with c = 1 in agreement with our expectation. We 
note that the intermediate plateaux to = 3/6 and to = 5/6 
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FIG. 7. (Color online) Block entropy S v n(£) vs block length d(£\L) 
(starting at one end of the tube) for several magnetization plateaux on 
L — 32 spin tube. Coupling constants are fixed to J\\/J± — 0.02. 




d(l\L) 



d(i\L) 



FIG. 8. (Color online) Block entropy 5V, at {€) vs block length d(£\L) 
(starting at one end of the tube) for L — 16 (open symbols) and 
L = 32 (filled symbols). Coupling constants are fixed to J\\/J± = 
0.1. (a) Magnetizations corresponding to plateaux; (b) intermediate 
magnetizations. 



also possess critical degrees of freedom, which could be com- 
patible with c = 1 or a slightly smaller value. We plan 
to investigate in the future whether there could exist a non- 
gaussian criticality nearby, or if it is simply due to numeri- 
cal uncertainty when coupling constants have very different 
amplitudes. Anyhow, for intermediate plateaux, there is no 
simple chirality language since more than two states per trian- 
gle (respectively four and three) are necessary to describe the 
low-energy configurations. 

Now, we can increase the coupling constant Jy / J± since 
the effective Hamiltonian and the CORE analysis indicate that 
the chirality degree of freedom should become gapped beyond 
some critical values (see Table I for S = 3/2). For instance, 
when fixing Jii /J± = 0.1, our data shown in Fig. 8 confirm 
that chirality has become gapped in the upper plateau, but re- 
mains gapless (with c = 0.97) for m — 1/6. The critical ratio 
that we find for gapless chirality on m = 7/6 plateau is close 
to 0.04, while chirality is always gapless on the m = 1/6 as 
long as it exists, i.e. for J\\/J± < 0.25. The quantum phase 
transition between gapless and gapped chirality phases for the 
upper plateau is indicated on the phase diagram (Figure 6). 
The central result is that both critical values are in excellent 
agreement with our range-2 CORE estimates (see Table I). 

For completeness, we also plot in Fig. 8(b) the scaling of 
the block entropy for intermediate magnetizations, which all 
correspond of course to critical gapless phases. In particular, 
at low magnetization, data are compatible with a 2-component 
Luttinger liquid with c = 2 as predicted. 25 For all the other 
magnetizations, our data are compatible with a single gapless 
mode c = 1. It is beyond the scope of this work to study the 
interplay between chirality and magnetic degrees of freedom 
outside magnetization plateaux, but it could be interesting to 
investigate the stability of the ferrochiral phase for arbitrary 
magnetic field. 45 



3. Nature of the gapped phase 

From the path-integral approach and the bosonization of 
the effective Hamiltonian, several different gapped phases are 
predicted to possibly occur when varying the coupling. How- 
ever, the block entanglement entropies does not give any in- 
formation about the nature of the gapped phase observed for 
J\\/J± > 0.04 in the upper m = 7/6 plateau. To investi- 
gate more precisely this question, we have computed the local 
magnetization values for different couplings, shown in Fig. 9. 

In a very narrow range close to the phase transition, a clear 
staggered spin imbalance is observed in our simulations (see 
Fig. 9a), as predicted from the bosonized Hamiltonian (27). 
The local magnetizations vary around their mean value m — 
7/6 with one chain having a clear different magnetization than 
the two others (in fact, we cannot exclude the possibility that 
all three chains will have different magnetizations). While 
the symmetry cannot be broken on a finite lattice, it turns out 
that DMRG simulations get locked in one of the degenerate 
ground-state. Note also that the level crossing found by ED in 
this region could impede the accuracy of DMRG results here. 
It is interesting to contrast our result with the small uniform 
spin imbalance phase found in the S = 1/2 tube which seems 
to signal the entrance into a regime where the pseudo-spin 1/2 
effective Hamiltonian is not valid anymore. 21 

When increasing slightly J\\/J±, but still deep in the 
plateau phase, we observe that all chains recover the same 
magnetization (see Fig. 9b), in agreement with having a fer- 
rochiral phase. 

In order to ascertain the validity of the description in terms 
of chiral degree of freedom, as described in Sec. II C, we have 
computed by ED on a small 3x6 cluster with PBC the weights 
of these degrees of freedom in the reduced density matrix of 
one triangle (see Refs. 39 and 40 for a discussion of this tech- 
nique). In a wide range of Jn / J± < 0.08 including the three 
different phases, XY, staggered spin-imbalance and ferrochi- 
ral, we find that the weights of these two states exceed 90%, so 
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FIG. 9. (Color online) Local magnetizations obtained by DMRG 
simulations on a 3 x 32 spin tube with spin 3/2 on the m = 7/6 
plateau. Upper and lower panels correspond respectively to (a) 
J||/Jl = 0.04 and (b) J n /J ± = 0.06. 



that we are rather confident that the effective model in terms 
of chirality remains valid. 

III. THE TWISTED SPIN TUBE 

A. The model and its experimental realization 

Experimentally, only a few materials have been suggested 
to realize spin tube geometries. One such geometry, which 
we will study in this section, corresponds to the com- 
pound [(CuCl2tachH)3Cl]Cl2. Magnetic measurements 46 
have shown that it forms a twisted triangular spin tube. We 
call it "twisted" because of the different structure compared 
to the simple tube of the Sec. II. The spins S = 1/2, coming 
from the copper ions are arranged in a one-dimensional array 
of equilateral triangles, and each spin of a triangle is coupled 
to the spins of the two others chains of the neighbouring tri- 
angles (Fig. 10). This corresponds to add diagonal couplings 
to the model (2) while the longitudinal one J\\ vanishes. The 
Hamiltonian describing the twisted spin tube reads 

H = H± + Hd + Hh 

j <*=1,2,3 



S a -i 



j+lj 



(28) 



H h 



j q=1,2,3 

-*E E s «,p 

j q=1,2,3 



where Jd is the diagonal antiferromagnetic coupling. This 
model is believed to describe [(CuCl2tachH)3Cl]Ci2 for the 
values J±_ = 0.9K and Jd = 1.95K 46 , thus the compound 
belongs neither the strong or weak coupling regime. Theor- 
ical and experimental investigations 4748 have shown that this 



compound behaves essentially as a spin-3/2 chain, displaying 
a Luttinger liquid behavior. 




FIG. 10. Lattice structure of the twisted spin tube. 



B. Path integral approach 

Following the same steps as in Sec. II B, we start by finding 
the classical ground state of the Hamiltonian (28). Classically, 
the extreme cases J± — and Jd = are easily understood 
in the absence of the magnetic field. For J± = 0, the lat- 
tice is bipartite and the ground state is the Neel state. On the 
other hand, for Jd = the triangles are decoupled and, as for 
the regular spin tube, the three spins are coplanar with a 27r/3 
angle between them. It turns out that the lowest-energy con- 
figuration is one of those two states in the whole range of the 
coupling parameters. 46 From Jd/J± = to Jd/J± = 3/2, 
the ground-state is the 2-7r/3 state and for higher values of 
Jd/J± it changes to the Neel state. This can be seen from 
a calculation in Fourier space by minimizing the resulting ex- 
change coupling J(kx,k\\) = 2cos(k±)(J± + 2J<jcos(fe||)), 
where fcj^ = 0, 27r/3. We have checked numerically that the 
effect of a magnetic field is just to polarize the spins in these 
two configurations and does not change the transition value 
Jd/J± = 3/2. Thus, to investigate this model using the path 
integral approach, we have to treat the two regimes separately. 



1. Regime J d /J±_ < 3/2 

In this range of Jd/J±, all the triangles are in the same 
state, with angles of 120° between neighbouring spins par- 
tially polarized by the magnetic field. This comes from the 
three-colorability of the lattice. This state is actually the same 
than the umbrella structure for the simple spin tube, but with- 
out the staggered order along the tube. We write this state as 



ij q _ j — iJ 



sm(6» )cos(v5°; 
sin(0 n )sin(y>°) 
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cos(6> ) 



where cos(6*o) 



h6X, 



and ifP a — (a ~ l)27r/3 up to 



S(3J ± 

an additional constant. Introducing the quantum fluctuations 
0a -> a ,j = ^o + S6 a ,j, <Pa -* fa + <Pa,j and the conjugate 
momentum H a ,j, we get the same expansion (8) for the spin 
operators, except for the alternate order term (— 1) J '. Then the 
action in the continuum limit reads 
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This action has the same form than (9) for the simple tube, 
the only difference appearing in boundary terms. We observe 
that, except for the longitudinal part in (d x ip a ) 2 , the couplings 
Jx and J d play the same role. More precisely, at large-scale, 
one can argue that the diagonal coupling J d is essentially iden- 



tical to the perpendicular one. The factor of two for J d simply 
tells that there is twice as many diagonal couplings than per- 
pendicular per unit cell. Thus we can follow the same steps 
and after performing the gaussian integration we find 



S[{fa v }] = S ch [fa,fa]+S a [fa] 

S*[<h,<h] = fdrdxi^X^ [{d T faf + (d T cf> 2 ) 2 ] + -A^ [{d x faf + (d x fa) 2 ] - ifi(fad T fa - fodrfa) (M) 

s s {cb s } = I d T dx^-\^\d T 4> s ) 2 + \\ x s Xd x <j> s ) 2 + id S ^d^ s Y 



where the constants are functions of the microscopic parame- 
ters. 

At this order the action is again decoupled into a symmet- 
ric action for <fi s and a chirality action for the two other fields 
(\>x and <j>2. Surprisingly, the mass of the latters vanishes af- 
ter momentum IT Q integration. However, it does not invali- 
date our previous discussion about the necessary rescaling of 
<j) s . Indeed, the expansion (8) for the spin operators and for 
the Hamiltonian is restricted to second order in the fields. Ex- 
panding up to fourth order, mass terms in (ip a — (fia+i) 4 would 
appear and we would recover a mass term, ensuring that fa 
and fa are still small. Appart from this mass cancellation, 
there is no difference with the simple spin tube case. This is 
expected as they both have the same classical configuration 
in this regime of J d /Jx- For the twisted tube, the staggered 
order is in the diagonal coupling, as the angle between spins 
coupled by J d is larger than ir/2. 

Then we perform the duality transformation on the sym- 
metric part of the action. We obtain the plateaux existence 
condition 3(5 — to) <E Z, which is obviously the same as for 
the simple tube since it does not depend on the detailed geom- 
etry but only on the unit cell. We also compute the expression 
of the cosine operator dimension 



\ 



3J d (S 2 - to 2 ) 



J\ 1 



2±l 



(32) 



which is very similar to the simple tube case (15) and bears 
the same functional form with in. From the dependence on the 
microscopic parameters, we predict again to observe plateau 
in the strong coupling regime along the rungs. 

Concerning the fields fa and fa, they describe the chirality 
degree of freedom as for the simple spin tube. The form of the 
action being the same, the same reasoning than in Sec. II B 3 
holds and the same four phases are possible. Particularly, we 
claim again that the Berry term in this action makes possi- 
ble a chiral gapless phase for some values of the microscopic 
parameters we are not able to compute (but still a priori in 
the strong coupling regime). Moreover, we expect that in the 
present case, there are more chances to be in this phase com- 
pared to the simple tube case because of the vanishing bare 



2. Regime J d /Jx > 3/2 

We start here from the partially polarized Neel state, 
parametrized as 



S, 



<*>j 



(-lpsin(6>o)cos(^ ) 

-■ : .S' I (-l)*sin(0o)sin(po) 

cos(6» ) 



(33) 
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where cos(#o) = h/(8SJ d ) and we choose ipo — 0, this free- 
dom of choice reflecting the U(l) degeneracy of the ground 
state. 

Then we proceed as usual, allowing these angles to fluc- 
tuate by small quantities 56 a j and (f a j as 9q — > a j — 



6q + S9 a j, <p Q a — > ip a .j, and introducing new variables H a .j 
conjugates of the V5 QJ 's. Following the same steps as pre- 
viously, namely rewriting the Hamiltonian as a function of 
these fluctuation variables, expanding up to second order in 
the fields and taking the continuum limit, we obtain the action 



s[{u a },{<p a }] = JdTdx( ]T 



aJ d {S 2 - m 2 )(d x ip a ) 2 + a(2J d - J ± ) 



S 2 



a=l,2,3 

1 , n r r ,S 2 -m 2 
+ -(2J d -J X ) 

z a 



S 2 



-IT 

,2 a. 



[(ip x - ip 2 ) 2 + (ip 2 - (f3) 2 + (f3 ~ ¥>l) S 



+ al2J d [ 1 

a=l,2,3 



S 2 — m 2 
S — m 



J±- 



S 2 



S 2 — m 2 
d T (p a - TL a d T (p a 



(34) 



(nin 2 + n 2 n 3 + n3n!) 



The condition 2 J d > J± for the action to be positive-definite 
is automatically fulfilled as we started from the assumption 
Jd/J± > 3/2. The colinear nature of the classical ground 
state is reflected in the absence of the terms IT Q (ip a +i—(p a -i) 



(see (9) or (30)). 

The next steps are to use again the transformation (11), per- 
form the gaussian integration in the II Q fields and rescale the 
symmetric field as <p s — > <^> s /v3. This leads to the action 



SK^a}] = S c h[(t>l,(t>2] + S s [(/) s } 

s[{<t> a }} - I ' dTdx\ l -\^ [(dr^) 2 + (d T fo) 2 } + h^ [(d^) 2 + (dxfo) 2 } + m 2 ^ 2 



S s [4> s ] = J drdxi l -\^{dAsf + \^ ] {d x 4> s ) 2 + i3^^d T <b s 



(35) 



Although we end with the same decoupling as previously, 
this action is actually quite different than (10) and (31). The 
crucial point is that there is no Berry term i(<pid T 4>2—(t>2dT<t>i) 
here. Thus those fields are automatically gapped, and there 
is no possibility of neither an emergent gapless phase nor 
a spin-imbalance phase (no z(II\I/* — n*^) term), contrary 
to the regime J d /J± < 3/2 or for the simple tube. This 
is not surprising however, given the fact that we are not in 
the strong-coupling regime and so there is no possibility for 
the chirality described by those two fields to be gapless (or 
more explicitly, there is no chirality in a colinear configura- 
tion). Then they only correspond to high-energy excitations 
and we can integrate them out by using the saddle point solu- 
tion (j>i = (p2 = 0. 



It remains only the action for the symmetric field, which 
is exactly the same we have already encoutered. We apply 
the duality transformation and repeat our analysis of the dual 
action. At the end, we recover the plateaux condition 3(S* — 
to) 6 Z. On the other hand, the situation is very different 
for the scaling dimension, as we find it to be independent of 
the parameters of the microscopic model in this second order 



calculation. It reads 



A=| 7 r V / 5^ 



(36) 



This means that a plateau will either be always present or al- 
ways absent when J d /J± > 3/2. More precisely, if a plateau 
is absent for a given spin value, then it will also be absent for 
higher spins. 



3. Discussion 

Using the results derived above, we are able to discuss 
the case of [(CuCl2tachH)3Cl]Cl2, which belongs to the 
Jd/J± > 3/2 regime. For S = 1/2, the OYA condition pre- 
dicts only one plateau at a magnetization per spin m = 1/6, 
or 1/3 of the saturation value. The scaling dimension is 
A = 7r/\/2 > 2 for these values, so the cosine is irrele- 
vant. With this result, we predict that the compound does not 
possess any plateau in its magnetization curve. In a previous 
work, Fouet et al. A1 reached the same conclusion. They found 
that for the realistic coupling values J d /J± = 2.16, the model 
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(28) behaves as an effective spin-3/2 antiferromagnetic chain 
with no plateau, and their DMRG calculations confirmed this 
effective Hamiltonian approach. 

Overall, our analysis is consistent with their numerical 
phase diagram based on the DMRG. Indeed, they observed a 
finite size for the m = 1/6 plateau from Jd/J± — to about 
Jd/J± = 3/2, where it vanishes. We find here the same re- 
sult. Altough the form of the action of the symmetric field is 
the same starting either from the 2tt/3 state for Jd/J± < 3/2 
or from the colinear Neel state for Jd/J± > 3/2, the scaling 
dimension of the cosine term indicates that in the first case 
the plateau is always present while it disappears in the sec- 
ond one. Thus, we expect qualitatively the same phase dia- 
gram with two different regimes. However, we should con- 
sider carefully this result, as for higher spin 5 = 3/2 DMRG 
calculations seems to indicate that the plateau m — 3/6 van- 
ishes also at Jd/J± = 3/2, or very close to this point (see 
below). Thus the important point is that we have obtained two 
different results for the cosine dimension, one being indepen- 
dent of the microscopic parameters. The other one tells us that 
we could observe plateaux depending on the value Jd/J±- 
But we should keep in mind that the critical values predicted 
come from a large 5 analysis. 



C. Strongly coupled chains 

The triangular unit cell being the same as in the simple tube 
case, the same procedure as done in Sec. II C can be applied 
to build an effective Hamiltonian on the lowest and the high- 
est magnetization plateaux. For instance, the effective Hamil- 
tonian to first order perturbation in Jd/J± for the spin- 1/2 
twisted spin tube on the unique in = 1/6 plateau reads 47 



H, 



eff 



T» 



•2(r. 



3 T H 



+ T i T U 



i)]> (37) 



where the chirality operator r is defined in (22). Thus, to go 
beyond this perturbation theory, we propose the same form 
(26) for general Hamiltonians describing the emerging chiral- 
ity degree of freedom along with its bosonized form (27). It 
follows that the same succesion of a gapless phase then two 
gapped phases is predicted, as for the simple tube. 

Keeping only the XXZ part of (26), we use again a range-2 
CORE calculation as in Sec. II C 2 to extract numerical values 
for J xy and J z . Surprisingly, at this level of approximation, 
we find a completely different result compared to the simple 
tube case, namely the absence of the XY-ferrochiral transition. 
For the twisted tube, the chirality remains in the | A | < 1 phase 
on both plateaux as the coupling J\\ increases. Moreover, we 
also find this absence of a gapped phase for higher half-integer 
spins, contrary of the simple tube study which shows the dis- 
appearance of the XY phase in the large-5 limit. From this 
CORE calculation we do no find a gapped phase. 



D. DMRG results for S = 3/2 

We now consider the 5 — 3/2 case using large-scale nu- 
merical simulations with DMRG algorithm (details are identi- 
cal to Sec. II D). First, we confirm that magnetization plateaux 
that satisfy OYA criterion 3(5 — m) G Z exist for small inter- 
triangle coupling, such as Jd/J± = 0.1 as shown in Fig. 11. 




FIG. 11. Magnetization curve (m) vs magnetic field h for the twisted 
three-leg spin tube in the case S = 3/2. DMRG simulations were 
performed with L — 32 and Jd/J± =0.1 

By performing a finite-size analysis of the plateaux widths, 
we can obtain the phase diagram in magnetic field, shown in 
Fig. 12. As for the simple tube, the precise location of the 
plateau regions can be quite hard due to BKT transitions. In 
the plot, we have used as a simple criterion that the extrapo- 
lated plateau should be larger than 0.005 J± to be considered 
as finite. Similarly to Fig. 6, there also exists a plateau at 
m = 0, which corresponds to a spontaneous dimerization of 
the tube and is not the subject of our present study. We find 
that the largest plateau corresponds to m = 1/3 and is stable 
in all the region Jd/J± < 3/2; it exhibits a small anomaly 
at the tip, similar to what is found in the 5 = 1/2 case. 47 
For a given Jd/J± and increasing h, we also note that the or- 
der of disappearance of the plateaux disagrees with the path- 
integral prediction, as in the simple tube case, which may be 
due to renormalization effects in the field-theory parameters 
since the field-theory is valid at large 5 and we are consider- 
ing 5 = 3/2 here. 

Following Sec. II D 2, we now turn to the investigation of 
the chirality for the extreme plateaux m = 1/6 and m = 7/6. 
According to the range-2 CORE Hamiltonian of the previous 
section, we expect to observe criticality for these degrees of 
freedom in all the plateau. Computing the scaling of the block 
entanglement entropy with DMRG (data not shown), we ob- 
serve that while it seems to be the case for m = 7/6, we do 
find a transition to a fully gapped regime in the lower plateau 
m = 1/6 when Jd/J± = 0. 12 (red hashed regions in Fig. 12). 
But it turns out that the mechanism gapping the chirality is 
here very different from the scenario determined for the sim- 
ple tube. Computing the local magnetizations for the three 
chains, we find neither a staggered spin-imbalance phase nor 
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FIG. 12. (Color online) Phase diagram of the twisted spin tube with 
S = 3/2 as a function of the couplings Jd and magnetic field h 
(«7j_ is fixed to 1). Several magnetization plateaux can be observed 
(filled regions) and an additional m = plateau is also found (but 
is much smaller and not shown on the figure). Data correspond to 
numerical simulations on 3 x 32 lattice with DMRG. Red hashed 
areas correspond to critical chirality regions (see text). 



a ferrochiral phase. Instead, we find for Jd/J± > 0.12 a 
unique uniform spin-imbalance following the critical phase. 
The imbalance is very strong as shown in Fig. 13, one chain 
having a negative magnetization, and it holds up until the dis- 
appearance of the plateau. 

The explanation for the absence of the two gapped phases 
present in the strong coupling limit is given by the reduced 
density matrix weights. In contrast to the simple tube case, 
the weights of the two chiral states rapidy becomes quite small 
(for instance 47% for a 3 x 6 twisted tube with Jd/J± = 0.3 
as found by ED). Thus, for such values we can not rely on the 
effective Hamiltonian. Note that a similar argument involving 
a mixing with other triangle states, thus prohibiting the use of 
the effective model, was given in Ref. 21 to explain the uni- 
form spin imbalance phase found for the simple S = 1/2 spin 
tube at magnetization m = 1/3. Instead we have to rely on 
the path-integral results, valid for any Jd/J±, which indeed 
predict this uniform spin-imbalance phase. The important re- 
mark here is that, despite the important difference between 
the two tubes, we are able to understand all the phases ob- 
served by combining the results of the path-integral and of the 
bosonized form of the strong-coupling Hamiltonian. 



IV. CONCLUSION 

In this paper we have studied the magnetic and non- 
magnetic properties of frustrated three-leg spin tubes under 
a magnetic field. We have considered two kinds of geome- 
tries, one of which is relevant for the recently studied com- 
pound [(CuCl2tachH)3Cl]Cl2. Our first result concerns the 
presence of plateaux in the magnetization curve. We give the 
value of the magnetization at which such plateaux can appear 




(a) 



(b) 



Site index i 

FIG. 13. (Color online) Local magnetizations obtained by DMRG 
simulations on a 3 x 32 twisted spin tube with spin 3/2 on the m = 
1/6 plateau. Upper and lower panels correspond respectively to (a) 
J d /J ± = 0.2 and (b) J d /J ± = 0.3. 



given the magnitude of the spins S, as well as the critical cou- 
pling for which such plateaux are expected to appear. We have 
used two complementary techniques. The first one is the path 
integral method which, because of the topological nature of 
the Berry phase term, gives trustable qualitative results for any 
value of S, and whose prediction of critical couplings are ex- 
pected to be also quantitatively accurate for large S. The sec- 
ond technique used is the DMRG method which is however 
more suited to relatively small spins S. 

While magnetization plateaux are not specific to frustrated 
systems, there are emergent low-energy degrees of freedom 
which presence is due to frustration. For historical reasons, 
we have dubbed those degrees of freedom chirality degrees of 
freedom. Their origin comes from degeneracies in the ground 
state for decoupled triangles, which motivates the use of a 
third technique to complement the path integral and DMRG. 
This strong coupling technique corresponds to studying an ef- 
fective Hamiltonian in a reduced Hilbert space were high en- 
ergy degrees of freedom are neglected. This supplementary 
degree of freedom can remain gapless even when the magne- 
tization degrees of freedom are gapped (in the magnetization 
plateau). Here again, the agreement between the three tech- 
niques for predicting the critical couplings at which such de- 
gree of freedom disappear is qualitatively excellent and quan- 
titatively quite satisfactory. It is important to stress that, al- 
though we restricted ourselves to the cases of two degenerate 
spin-1/2 representations for half-integer spin cases, the chiral- 
ity is in principle a generic feature that can give rise to more 
complicated effective Hamiltonians. It could arise also for in- 
teger spin tubes, provided they are tuned to the appropriate 
value of the magnetization. 

The results obtained here and the excellent complementar- 
ity of the path-integral technique, the effective hamiltonian 
approach and the DMRG calculations are very encouraging. 
Indeed, the study of potentially gapless non-magnetic degrees 
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of freedom has become a central topic in the study of exotic 
phases in frustrated quantum magnetism. In the systems an- 
alyzed here, one important generalization that deserves fu- 
tures studies is the interplay of doping with such degrees of 
freedom. It is by now relatively well established that doping 
will result in shifts and splitting of the magnetization plateaux 
in such quasi one-dimensional systems 49-51 , but how doping 
may affect the non-magnetic degrees of freedom is still an 
open problem. The other important extension concerns the 
role of such non-magnetic degrees of freedom in higher di- 
mensional frustrated systems, where some results for distorted 
kagome lattices are indeed quite instructive. 52 ' 53 A chirality 
also appears for instance in the study of trimerized Mott insu- 
lators. 54 The emergence of new analytical and numerical tech- 
niques for studying such issue in two and three dimensional 
frustrated magnets are also very promising. 
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where we have defined J = (J T + 3 s ~ m , J x ). The vorticity- 
free part simply leads to a zero divergence constraint on the 
auxiliary field d T J T + d x J x — and we obtain 



S[<ps,v, J] = / drdx 



I 2Ai s) V 



-3- 



.S- 



H tztJx + iJ^d^s.v 



2A 



« 



(A3) 



The constraint can be solved in one dimension by introducing 
the dual field $ s defined by J M = t^d v ^ s , this field beeing 
vorticity-free. Then we integrate by parts the last term in (A3) 
and, with the redefinition $ s = $ s — 3 s ~ m x, we get 
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Appendix A: Duality transformation 



S[$ a ] = f 'dTd&j— ^-(0 r * 



i2irp v $ 



2\ { T s) 



(d x $ s 



S-m 
6 x 



(A4) 



We present here the details of the duality transformation 
used on the symmetric part S s of the action (10). First, we 
perform a Hubbard-Stratonovich transformation introducing 
an auxiliary field J = (J T ,J X ) and we divide the field <f> s 
in two parts <f> 8 — <j> SyV + <j) s j, where <j> s j has no vorticity 
{d^d u — dyd^)4> s j = 0. The action reads 



S[4> s ,v,4>sj,J] = / drdx 



1 



2\ { T s) 



31 



1 



2A£ S) 



-/ 



i [ J T + 3 ) d T <fi s . v + iJ x d x (j> s . v 

i [ J T + 3 ) d T <l) s j + iJ x d x <p, f 



(Al) 

Integrating by parts the last two terms containing the vorticy- 
free component, the action takes the form 



s[4>s,v,4> s j,J} 



drdxl — rT J T 
l2Ai s) V 



S — m 



1 



+ — j-cJZ + i(J T d T + J x d x 



i(d T J T + d x J : 



x)<Ps,f 



In this action p v is the space-time density of vortices de- 
fined as (d T d x - d x d T )4> v = e^d^d^ = 2kJ2j 1j,-uK r ~ 
Tj lV )S(x — Xj, v ) — 2irp v with (tj v ,Xj^ v ) the space-time coor- 
dinates of the j-th vortex and qj, v € Z its charge. After sum- 
ming over all the vortex configurations in the partition func- 
tion and rescaling the imaginary time, we finally end with the 
action 



S[*.] = / drdx{ ^(Vl> s ) 2 



gicos 2tt 



i n S — m 
f. + 3 x 



(A5) 



where K = 1/ 

we have not calculated. 



\ { r ] \ { i\ V 



(d T , d x ) and g\ is a constant 
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